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This paper reports ah-atom computer simulations of five phospholipid membranes (DMPC, 
DPPC, DMPG, DMPS, and DMPSH) with focus on the thermal equilibrium fluctuations of volume, 
energy, area, thickness, and chain order. At constant temperature and pressure, volume and energy 
exhibit strong correlations of their slow fluctuations (defined by averaging over 0.5 nanosecond). 
These quantities, on the other hand, do not correlate significantly with area, thickness, or chain 
order. The correlations are mainly reported for the fluid phase, but we also give some results for the 
ordered (gel) phase of two membranes, showing a similar picture. The cause of the observed strong 
correlations is identified by splitting volume and energy into contributions from tails, heads, and 
water, and showing that the slow volume-energy fluctuations derive from van der Waals interactions 
of the tail region; they are thus analogous to the similar strong correlations recently observed in 
computer simulations of the Lennard- Jones and other simple van der Waals type liquids [U. R. 
Pedersen et al., Phys. Rev. Lett. 2008, 100, 015701]. The strong correlations reported here confirm 
one crucial assumption of a recent theory for nerve signal propagation proposed by Heimburg and 
Jackson [T. Heimburg and A. D. Jackson, Proc. Natl. Acad. Sci. 2005, 102, 9790-9795]. 



In 2005 Heimburg and Jackson showed that biomem- 
branes may carry solitonic sound waves whose maximum 
amplitude and minimum velocity are close to the propa- 
gation velocity in myelinated nerves . Their paper con- 
cluded: "It would be surprising if nature did not exploit 
these features." Subsequent works by the same authors 
argue directly that nerve signals are not primarily elec- 
trical, but solitonic sound waves carried by the nerve cell 
membrane O |3] . The conventional wisdom is that nerve 
signals propagate via electrical current as formulated in 
the Hodgkin-Huxley theory [4]. The Heimburg- Jackson 
theory explains anaesthesia as a straightforward effect of 
melting-point depression of the order-disorder transition. 
This can in turn explain why biomembranes have a tran- 
sition temperatures close to physiological temperatures 
|5i. 

The Heimburg-Jackson nerve-signal theory motivated 
this study. One element of the theory is the assump- 
tion that volume and enthalpy correlate strongly in their 
thermal equilibrium fluctuations, at least as it regards 
the slow parts of these fluctuations. This assumption 
was justified by experimental observations that mem- 
branes specific heat and compressibility are proportional 
in their temperature dependence across the phase tran- 
sition [6| |7]- Also, a correlation between lipid area and 
enthalpy has been proposed [6 based on experimental 



observations [8^. If so, both volume, membrane area and 
enthalpy are controlled by a single parameter. Below, we 
present first results from extensive computer simulations 
of different phospholipid membranes performed in order 
to investigate whether strong volume-energy and/or area- 
energy correlations are observed and, if so, what causes 
them. 

Recently, we studied equilibrium thermodynamic fluc- 
tuations of much simpler systems, namely various model 
liquids like the standard Lennard- Jones (LJ) liquid and 
similar systems [H [TOl [TTJ [121 • many such simple 
liquids one finds a strong correlation between the equi- 
librium fluctuations of virial W and potential energy [/, 
when fluctuations are studied at constant particle num- 
ber N ^ constant volume V ^ and constant temperature 
T (the so-called NVT ensemble [M]). Recall [15], that 
the virial W = W{t) gives the non-ideal contribution 
to the instantaneous pressure p = p{t) via the defining 
equation: p{t)V = NkBT{t) + W{t) where T{t) is the 
instantaneous temperature defined via the instantaneous 
kinetic energy. The virial is a function of the particle po- 
sitions [15 . For the LJ liquid, as well as for a united-atom 
toluene model, a dumbbell model, the Kob- Andersen bi- 
nary LJ liquid [16 , and other van der Waals liquids, W 
and U correlate better than 90% in their equilibrium fluc- 
tuations [9 . This reflects an effective inverse power-law 
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potential dominating fluctuations, as discussed in details 
in Refs. |TTJ [121 Liquids with poor W — U correlation 
include water and methanol [H [11]. In these cases the 
correlations are ruined by the hydrogen bonds, which are 
conventionally modelled via Coulomb forces - the exis- 
tence of competing interactions prevents strong W — U 
correlation in hydrogen-bonded liquids. 

For liquids with time-scale separation, like highly vis- 
cous liquids, strong W — U correlations are particularly 
signiflcant: Viscous liquids with strong W — U correla- 
tions are close to being "single-order-parameter liquids" 
in the classical Prigogine-Defay [13 [TBI sense [TOl [T51 [T9] . 
This implies that complex frequency-dependent thermo- 
viscoelastic response functions like the isobaric/isochoric 
dynamic speciflc heat, the dynamic thermal expansion 
coeflicient, the dynamic compressibility, etc, are all given 
by a single function [19]. In particular, these cannot "de- 
couple" [20] from one another - they must all exhibit 
relaxations in the same frequency range. It has also been 
shown that strongly correlating viscous liquids obey den- 
sity scaling, i.e., that if the relaxation time r is measured 
at different temperatures and density p, r is a unique 
function ^ of p^/T where the exponent 7 may be de- 
termined from studies of W — U correlations at a sin- 
gle state point [22]. Finally, it was recently found that 
strongly correlating viscous liquids have much simpler 
aging properties than viscous liquids in general [23] . 

Fluctuations are ensemble dependent, of course, and 
one may ask what happens if fluctuations are studied in- 
stead in the ensemble of constant temperature and pres- 
sure {NpT ensemble). In this ensemble virial fluctuations 
are not interesting, but there are strong correlations for 
simple liquids between the fluctuations of volume and po- 
tential energy [10]. This is the ensemble used below for 
studying biomembrane thermodynamic fluctuations. 

Motivated by the findings for simple model liquids and 
the Heimburg- Jackson theory, we decided to focus on 
phospholipid membranes in order to investigate whether 
the Heimburg- Jackson assumption of strong volume- 
enthalpy correlations is confirmed. A phospholipid has 
van der Waals interactions between its acyl chains and 
hydrogen bonds in the head region. Similarities to sim- 
ple van der Waals liquids are not at all obvious, and 
the microscopic origin of the volume-enthalpy correla- 
tion tentatively derived from experiments [6] [7] is not 
trivial. Thermodynamic fiuctuations of simulated mem- 
branes have been studied in the past [24l [25l |26l [27) . 
however, not with a focus on volume-energy correlations. 

Because of the large amount of water in the sys- 
tem and the hydrophilic head groups, one does not ex- 
pect strong instantaneous correlations of phospholipid 
membrane thermodynamic fiuctuations [9] [TTj. The 
Heimburg- Jackson theory, however, relates to strong cor- 
relations of slow degrees of freedom of the biomem- 
brane (on millisecond time scales), and so do the exper- 
iments they quote indicating strong correlations [6l [7]- 




FIG. 1: Snapshots of DMPC membranes in the fluid phase 
(DMPC-f; panel A) and in the ordered phase (DMPC-g; panel 
B). The red atoms are the oxygen atoms of water molecules; 
hydrogen atoms were removed for visual clarity (but included 
in the simulations). Acyl chains are colored green. The frame 
indicates the periodic boundary box. 

This is analogous to the situation for highly viscous liq- 
uids where time-scale separation between the fast, vi- 
brational degrees of freedom and the much slower con- 
figurational ones is also e.g. essential for viscous fiow 
or visco-elastic responses. Phospholipids are the major 
constituent of biological membranes [28 . Close to phys- 
iological temperature these membranes undergo a tran- 
sition from an ordered phase (L^) to a disordered phase 
(La) [3 [H [291 [301 [31]. Below, we evaluate the strength 
of V — U correlations for both the disordered and or- 
dered phases with main focus on the disordered phase. 
The correlation strengths are calculated for a range of 
time scales. We show that V — U fiuctuations correlate 
strongly, but only on long time scales [32]. We also in- 
vestigate how well membrane area as well as chain order- 
parameter fiuctuations correlate with V and U fiuctua- 
tions; such correlations are generally weak. Finally, the 
cause of the correlations is identified by splitting vol- 
ume and energy into contributions from the tail-, the 
head-, and the water region: The slow, strongly corre- 
lating V — U fiuctuations are shown to derive from the 
tails which are dominated by van der Waals interactions, 
thus establishing a conceptual link to the strong cor- 
relations of the slow pressure-energy fiuctuations - the 
virial/potential energy correlations - of simple van der 
Waals liquids [9]. 

SIMULATION DETAILS 

Details of the seven membrane simulations performed 
are listed in Table |T| The following abbreviations are 
used: DMPC-f and DMPC-g: fully hydrated di- 
myristoyl-phosphatidyl-choline membrane in the fiuid 
and ordered phases, respectively [33]. Figure [T|V show 
a configuration of the DMPC-f membrane. DPPC-f 
and DPPC-g: fully hydrated di-palmitoyl-phosphatidyl- 
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TABLE I: Overview of simulation details and results 



(A) System (phase) 


^sim N t^rod %pid ^ [K] ^actual 1^] ^waterMipid 


DMPC-f ("flnirl^ 
j_vivijr v_y i 1 iiU-iu. ) 




151 


121 


128 330.0 329.0±1.6 


33 


DMPC-g (ordered) 




65 


36 


64 286.0 285.3±2.0 


33 


DPPC-f (fluid) 




180 


124 


72 325.0 324.0±2.1 


29 


DPPC-g (ordered) 




78 


48 


64 304.0 303.2zb2.1 


33 


DMPG (fluid) 




149 


49 


128 330.0 329.0±1.6 


33 


DMPS (fluid) 




139 


49 


128 340.0 339.0±1.7 


36 


DMPSH (fluid) 




136 


35 


128 340.0 339.1±1.6 


37 


(B) System (phase) 




7 


Ruut 


RuA RvA RaSoo ^USon 


^VSnn 


DMPC-f (fluid) 




n 77-1-0 DzL 


c/.U 


0.82 


0.50 0.57 -0.75 -0.49 




DMPC-g (ordered) 


0.47=b0.07 


6.1 


0.31 


0.02 0.05 -0.64 0.12 


0.14 


DPPC-f (fluid) 




0.87=b0.04 


10.2 


0.89 


-0.29 -0.36 0.00 -0.61 


-0.71 


DPPC-g (ordered) 


0.75±0.07 


6.6 


0.71 


-0.16 0.12 -0.67 0.09 


-0.07 


DMPG (fluid) 




0.82±0.07 


8.5 


0.80 


0.41 0.40 -0.76 0.01 


0.08 


DMPS (fluid) 




0.59=b0.07 


7.6 


0.64 


0.30 0.28 -0.71 0.04 


0.20 


DMPSH (fluid) 




0.78±0.07 


13.2 


0.84 


0.43 0.51 -0.50 0.05 


0.14 



(A) tgjj^: Total simulation time in nanoseconds, t-^-^^^: Length of production run in nanoseconds (only membranes in 
quasi-equilibrium, i.e., with no detectable drift in the area per molecule, were included in the data analysis); A^ijpj^j: Number 

of lipid molecules; T: Thermostat temperature in Kelvin; ^actual- Average temperature and standart diviation of 
fluctuations; A^water/^lipid' Number of water molecules per lipid molecule. (B) Rfjv'. Energy- Volume correlation coeflicient 
(see Equation [T]) . The bar indicated a 0.5 ns average. The uncertainty is estimated from the DMPC-f and DPPC-f 
trajectories as described in the text (67% confldence interval); 7: Energy- volume scaling factor in units of 10~^ mL/J (see 



Equation 11); Ruijt'- Energy- "Energy of tail region" correlation coeflicient. Rua- Energy- Area correlation coeflicient; Ri^ 



Volume- Area correlation coeflicient. Rascd ' Area- "chain order-parameter" correlation coeflicient; i?^^^^ : Energy- "chain 
order-parameter" correlation coeflicient; Rysco ' Volume- "chain order-parameter" correlation coeflicient. 



choline membrane in the fluid and ordered phases, re- 
spectively [34 • DPPG: a fully hydrated di-palmitoyl- 
phosphatidyl-glycerol membrane in the fluid phase with 
calcium counter ions [35 . DPPS: a fully hydrated 
di-myristoyl-phosphatidyl-serine membrane in the fluid 
phase with calcium counter ions [35]. DMPSH: a fully hy- 
drated and protonated di-myristoyl-phosphatidyl-serine 
membrane in the fluid phase [35 . Phospholipids 
(and counter ions) were modeled using the all-atom 
CHARMM27 force-field [36^ with modified charges of the 
head group as described in detail elsewhere [33 . Mem- 
branes were hydrated using explicit water represented by 
the TIP3 model [37 . 

Simulations were carried out using the NAMD soft- 
ware package [38]. In all simulations a time step of 1.0 
fs was used. Temperature and pressure were controlled 
by a Langevin thermostat (damping coefficient: 5 ps~^) 
and a Nose-Hoover Langevin barostat (anisotropic reg- 
ulation; piston oscillation time: 100 fs; damping time: 
50 fs) [14. Electrostatic interactions were evaluated us- 
ing the Particle-Mesh-Ewald method [39] [40] with grid 
spacing of about lA updated every 4 fs. Lennard- Jones 
potentials were brought smoothly to zero by multiplying 
a switching function S{rij); S{rij) = 1 for Vij < i^on, 
S{nj) = {RlQ-rl)\Rl^^2rl-?>Rl^)/{Rl^-Rl^f 
for Ron < Uj < R^q and S{rij) = for > R^q 



with Ron = 10 A and RqQ = 12 A. Periodic boundary 
conditions were applied in all three dimensions. 

Initial configurations of DMPC-f, DPPC-f, DPPG, 
DPPS, DMPS and DMPSH are taken from Refs. [33l 
134) [35I Initial configurations of the ordered membranes, 
DMPC-g and DPPC-g, were build from a membrane sim- 
ulated by Venable and co-workers [41]. Throughout the 
simulation the acyl chains remain in an ordered structure 
as shown in Figure [l^ [31 . Thermal equilibrium was en- 
sured by monitoring the membrane area; only trajectories 
with no drift (compared to the thermal fluctuations) were 
used in the data analysis. Lengths of equilibrium trajec- 
tories are listed in Table jijA in the column under ^p^-od* 
The importance of firmly establishing that the systems is 
in equilibrium should be emphasized, since an apparent 
strong correlation would appear if volume and energy re- 
lax from some (arbitrary) out-of-equilibrium state. Only 
small finite size eflPects are expected at the simulated sys- 
tem sizes [42]. 

RESULTS 

The following collective quantities were evaluated ev- 
ery 0.5 ps: Potential energy /7, simulation box volume 
V = XYZ (where X, Y and Z are the box dimensions), 
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projected membrane area A = XY, box thickness Z, 
and average chain order parameter {ScD)ch- The lat- 
ter quantity characterizes the overah order of acyl chains 
[43] and is defined by {ScD)ch = 1(1 cos^ (Ocd) 
where Ocd is the angle between the membrane normal 
(z) and the C-H bond of a given methylene group, and 
{. . .)ch denotes an average over all methylene groups in 
all chains. 

First, we consider the fiuid DMPC-f membrane shown 
in Figure [T]A.. Instantaneous fiuctuations of volume and 
energy do not show any significant correlation (data not 
shown). This is not surprising, since a significant part of 
the simulation box is water, and water is known not to 
correlate [9 . However, if fiuctuations are averaged over 
time windows of 0.5 nanosecond, volume and energy cor- 
relate strongly as shown in Figure [2]A.. This is quantified 
by the correlation coefficient. 



{AUAV) 



uv 



V((Af/)2){(AF)2 



0.77, 



(1) 



where the bar here and henceforth indicates a 0.5 ns av- 
erage. R = corresponds to no correlation, whereas 
1 corresponds to perfect correlation. We will refer 



\R\ 



to R > 0.75 as strong correlation. For comparison, the 
correlation coefficient without averaging (Ruv) is 0.35. 

It is an appealing idea to establish a direct connec- 
tion between thermodynamic properties and microscopic 
structures, and several studies have focused on the or- 
dering of the acyl chains as the important microscopic 
structure |^. Following this philosophy, one pos- 
sible explanation for the observed strong correlation is 
that the order of the acyl chains is the single parameter 
controlling the ffuct nations: If chains, as a result of a 
thermal ffuctuation, become more ordered, one expects: 
decrease of energy, /7, volume, and area, A, but in- 
crease of thickness, Z, and {ScD)ch- We find R^a = 0-50 
and RyA = 0.57. Thus, the correlation has the expected 
sign, but it is significantly weaker than Rijy. The same is 
the case for Z and {ScD)ch' Figure [2^ shows the full cor- 
relation matrix. Clearly, a single parameter description 
is not sufficient. We need two parameters to describe the 
ffuctuations: one parameter controlling V and U and one 
"geometrical" parameter controlling A, Z, and {ScD)ch- 

Table ^ shows that DMPC-f, DPPC-g, DPPC-g, 
DMPG, DMPS, and DMPSH ah have strong volume- 
energy correlation {Ruy > 0.75). The volumes of 
DMPC-g and DMPS show weaker correlation with en- 
ergy. To show that this finding is not due to uncertainty 
from random noise, we estimated the error bar as fol- 
lows: Two of the simulations, DMPC-f and DPPC-f, are 
about three times longer than the remaining five simu- 
lations. We used these two long simulations to estimate 
the error bar of the calculated correlation coefficients by 
dividing them into three blocks of 40 ns (regarded as un- 
correlated blocks) ^44|. The standard deviation of the 
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1 : Energy 



2: Volume (XYZ) 



3: Area (XY) 



4: Thickness (Z) 



5: (Sr.r) 




FIG. 2: Correlations in the slow thermal equilibrium fluctu- 
ations of volume and energy (top) and correlation matrix for 
the DMPC-f membrane (bottom). The normalized fluctua- 
tions of volume and potential energy shown are averaged over 
time intervals of 0.5 nanosecond. Data are shifted and scaled 
such that the average value is zero and the standard deviation 
is unity. A significant correlation is observed and quantified 
by the correlation coeflficient, Rov — 0-77. This strong cor- 
relation can be associated with the tail region as seen by the 
similarities with Figure [5] The bottom panel represents the 
absolute values of the elements of the correlation matrix of 
energy, volume, membrane area, thickness and the average 
chain order-parameter where dark red illustrates strong cor- 
relation. Membrane area, thickness and average chain order- 
parameter are strongly correlated, but these quantities only 
correlate weakly with energy and volume. Similar results are 
found for the other fiuid membranes (except for DMPC-g and 
DMPS where the energy- volume correlation is only 0.47 and 
0.59, respectively), see Table |l] 



correlation coefficient is 0.03 and 0.10 for DMPC-f and 
DPPC-f, respectively - or about 0.07 on average. Each 
of the five short runs corresponds to a single block, and 
thus we estimate the error bar of these to be about 0.07 
(within a 67% confidence interval). For the two longer 
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runs we use the rule and estimate an error bar of 
0.07/y^3)=0.04. Thus, the uncertainty is smaher than 
the spread amongst i^^/y's, and the weak V-U correla- 
tion of the DMPC-g and DMPS membranes is genuine - 
a point returned to later. 

The term "slow fluctuations" has so far been defined 
via the 0.5 ns averaging time window (indicated with 
a bar). A more general approach is to investigate the 
following time-dependent correlation coefficient 



{AU{0)AV{t)) 



^{AU{0)AU{t)){AV{0)AV{t)) 



(2) 



Similarly, one defines the time-dependent energy-area 
correlation coefficient TuAit)- Figure [s] shows Tuvit) 
and TijA{t) for all seven investigated system. In contrast 
to TuA{t)^ we observe strong correlation on long time 
scales for TuY{t) {TuY{t) > 0.75) for the ffve membranes 
where slow ffuctuations of V and U correlate strongly 
consistent with R^y > 0.75 (Table |l^). In the following 
section, the slow parts of the volume and energy ffuctu- 
ations are investigated via the autocorrelation functions 
(AV(O)AV(t)) and {AU{0)AU{t)). 



Locating the slow volume and energy fluctuations 
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FIG. 3: Time-dependent correlation coefficient of volume- 
energy Tuvit) (left) and area-energy TuAit) (right). The 
fast fluctuations of volume and energy, t < 10 ~^ ns, only 
correlate weakly, whereas the slow fluctuations, t 1 ns, cor- 
relate strongly. Area and energy fluctuations are only weakly 
correlated in the investigated time regime. 



A membrane is a highly heterogeneous system , and 
it is reasonable to divide it into regions [27 . In the follow- 
ing three regions are deffned: t, /i, and where t (tail) 
refers to the hydrophobic acyl-chain atoms (i.e., atoms 
of methylene- and methyl groups in the acyl chain), h 
(head) refers to the hydrophilic lipid atoms (the remain- 
ing of the lipid atoms), and w refers to the water atoms 
(and counter ions). To identify the origin of the slow vol- 
ume ffuctuations, we construct Voronoi polyhedra [45] 
of heavy atoms (i.e., ignoring hydrogen) and sum the 
Voronoi volumes for the regions t, h and w. In this way, 
the total volume of the simulation box is divided into 
three terms, 

V = Vt^Vh^V^. (3) 

The auto-correlation function of the volume in Equation 
[2] can now be split into a sum of three auto- and three 
cross-correlation functions, 

{AV{0)AV{t)) = {AVt{0)AVt{t)) 

+ {AVh{0)AVh{t)) 

+ {AV^{0)AV^{t)) 

+ 2{AVt{0)AVh{t)) 

+ 2(AT4(0)AK,(t)) 

+ 2{AVh{0)AV^{t)). (4) 

Figure |4|\ shows these six functions for the DMPC- 
f membrane. The only nonvanishing function at long 



times (responsible for the slow ffuctuations) is the auto- 
correlation function of the hydrophobic (tail) part of 
the membrane, {AVt{0)AVt{t)). This is quantiffed by 
RvVt =0-94 being close to unity. 

The short time cross correlation between water and 
the tail region is signiffcant and negative. This is, how- 
ever, probably a spurious consequence of the Voronoi con- 
struction between neighboring groups of different size: 
When a water molecule approaches the boundary of the 
tail region it gains an (unfair) increase of the Voronoi 
volume, due to methylene groups being larger than wa- 
ter molecules. This show up as a negative correlation 
since the water molecule "steals" volume from the methyl 
groups, however, only at short times {t < 0.1 ns). 

In the simulation, the potential energy of the system 
consists of a sum of Lennard- Jones terms. Coulomb pair 
energy terms, and intramolecular binding-energies: 

where l/Jj"^ = Aeij{{(Jij /vijY'^ — {cFij/rij)^) are Lennard- 

Jones terms (van der Waals interactions), Uf^^^^ = 

QiQj / {Aireorij) are Columbic terms and U^^^^^ is a sum of 
the reminder terms (which are intramolecular; harmonic 
bonds, angles, dihedrals and improper dihedrals). Again, 
we split the total potential energy into contributions from 
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FIG. 4: Time- dependent auto- and cross-correlation function of volume and energy terms of the DMPC-f membrane. Panel 
A: Volume correlations functions of the tail-, head- and water regions (Equation |4|. Slow volume fluctuations are dominated 
by the tail region. The shot-time tail- water anticorrelation is probably an artifact of the way the Voronoi polyhedron is 
constructed (related to the tail- water interface, see text for details). Panel B: Correlation functions of energies of the three 
regions (Equation |8| . Again, slow fluctuations are dominated by the tail region. Panel C: Energy correlations split into 
intramolecular, Lennard- Jones (LJ) and Coulombic energies (Equation [7^ with x =tail). 



regions of tails, heads and water, 

U = Ut^Uh^U^ (6) 

where 

i = xj = ah i = xj= ah 

(7) 

with X either tails, heads or water. 

As for the volume, the auto-correlation function of the 
energy fluctuations in Equation [2] is also split into a sum 
of three auto- and three cross correlation functions, 

(AC/(0)AC/(t)) = (A/7i(0)Af/i(t)) 

+ {^Um^Uh{t)) 

+ {AC/™(0)A[/™(t)) 

+ 2{AUt{Q)AUh{t)) 

+ 2(Af/((0)AC/^(i)) 

+ 2(AC/^(0)AC/»(t)). (8) 



Figure [4^ shows the six auto- and cross correlation 
functions of /7t, Uh and U^. Again, the slow fluctuations 
are dominated by the tail region. It should be noted, 
though, that the slow tail-tail correlation is not quite as 
dominating, as for the volume fluctuations. This is quan- 
tified by Rijij^ = 0.82 (|l^) not being as close to unity as 
Ryy^ = 0.94: thus, slow energy fluctuations of the head- 
Uh and water region Uw are significant (in contrast to 
slow volume fluctuations of the head and water region). 

Figure |4p shows the auto- and cross-correlation func- 
tions corresponding to an additional splitting of the tail 
energy into "intramolecular interactions", "Coulombic 
interactions" and "van der Waals interactions" given in 
Equation [?) The van der Waals energies dominate the 
energy fluctuations of the tail region. 



DISCUSSION 

In the previous section, it was demonstrated that slow 
energy-volume correlation (of some membranes) origi- 
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FIG. 5: Normalized fluctuations of Voronoi volume and van 
der Waals energy of the tail region of the DMPC-f membrane; 
Vt{t) and Ui"^{t). Vt{t) and Ut^-^it) are shifted and scaled so 
the average value is zero and the standard deviation is unity. 
The correlation is strong with a correlation coefficient of 0.87 
(justifying Equation |9|. Note the similarity with Figure [2j\, 
thus slow V-U fluctuations are dominated by the tail region. 



nates from van der Waals interactions in the tail region; 
thus: 



AVf cx AU, 



LJ 



where 



LJ 



i=tailsj=all 



LJ 



(9) 



(10) 



This is consistent with our findings for simple strongly 
correlating liquids [Tl] [12] . It should be remembered, 
however, that Equation [9] neglected the energy terms Uh 
and Ut with some significance. If we correlate the slow 
fluctuations (defined as a 0.5 ns average) of the Voronoi 
volume of the tail region and the van der Waals energy 
of the tail region, we find a correlation coefficient of 
Ry^ijLj =0.87, as shown in Figure [5j The same num- 
ber for the fluctuations for whole simulation box was 
Ryij =0.77. The loss of correlation is associated with 
the neglected energy terms. Consistent with this. Table 
[ip shows a general correlation between Ryij and Ruu^^ 
e.g. for the weakly correlated membranes, DMPS and 
DMPC-g, both Ryij and RfjUt have low values. Appar- 
ently, slow energy fluctuations of the head region ruins 
the V-U correlation. A better understanding of how this 
loss of correlation is related to various lipid head groups 
deserves future study. 



It is convenient to define a volume-energy scaling factor 
via 



((AC7)2) 



(11) 



SO that if R 



uv 



1 then AV = jAU. For the seven 
membranes, given in Table [ij we find that 7 ranges 
from 6.1 X lO""^ mL/J (DMPC-g) to 13.2 x lO""^ mL/J 
(DMPSH). This scaling factor is indirectly accessible 
though experiments: Close to the phase transition tem- 
perature curves of excess enthalpy AH{T) and volume 
AV{T) can be superimposed [7 . Heimburg showed that 
such a relation is trivial if excess volume and enthalpy 
of relevant micro states are proportional [6 : AV = 
^ynAH. The inverse V-H scaling factor can be written 
as {^vh)~^ = tf = Since p = 1 atm= 10"^ 

J/mL is much smaller than 10^ J/ml (a typical 

value as seen in Table [l| then ^ and the pV con- 

tribution to enthalpy is negligible. Moreover, the excess 
volume and enthalpy of the phase transition is expected 
to be related to the slow equilibrium fluctuations close to 
the transition. Thus, the scaling factor calculated from 



equilibrium simulations (Equation 11) is closly related 



to the (experimental) "phase transition" scaling factor 
^yH and we expect them to have the same value. In 
agreement with this, Ebel, Grabitz and Heimburg re- 
ports a scaling factor of (9.01 ± 0.71) x 10~^ mL/J and 
(8.14 ± 0.67) X 10~^ mL/J for large unilamellar vesicles 
of DMPC and DPPC, respectively - only slightly smaller 
than our findings of 9.6 x 10~^ mL/J and 10.2 x 10~^ 
mL/J, respectively (Table |l^). This supports that we 
have identified the microscopic origin of the scaling of 
the phase transition. 

We find strong energy- volume correlations in both the 
fluid and ordered phases. The proportionality constant 
7, however, is not the same in both phases (Table [l^), 
but is about 2/3 of the fluid value in the ordered phase. 
Ebel, Grabitz and Heimburg [7 reported a similar de- 
crease (65%; multilamellar vesicles of DMPC or DPPC) 
when comparing scaling factors of the pretransition with 
the main transition. This suggests that the configura- 
tion of the ordered membranes (DMPC-g and DPPC-g) 
is related to the pretransition, and that the nature of the 
correlations changes when passing the phase transition. 
Interestingly, this is in contrast to our earlier findings 
for crystallization of the standard Lennard- Jones liquid 
[111 [12] . More work is needed to clarify the cause of this 
difference between biomembranes and simple liquids. 

Biologically relevant membranes are mixtures of many 
sorts of amphiphilic molecules and are far more com- 
plex than the membranes studied here. The origin of 
the strong correlations is, however, not specific, but re- 
flects the van der Waals bonded nature of the core of 
the membranes. Thus, biologically relevant membranes 
most likely also exhibit strong correlations. How solutes 
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in membranes effects the volume-energy correlation is an 
interesting subject for a future study. Another concern 
could be the finite size of the simulated membranes, since 
fluctuations on length scales larger than the length of 
simulations box cannot be represented ^26i. Again, since 
the origin of strong V-U correlations is intrinsic to van 
der Waals interactions [9l [10] , the correlations is not ex- 
pected to be affected by the size. 

One caveat is that although we do observe some cor- 
relation between lipid area and energy fluctuations, the 
correlation is weaker than the volume-energy correlation. 
This does not agree with the conjecture of Heimburg 
and coworkers [H E [6]. Their conjecture depends on 
the membrane being close to the phase transition, and, 
it is possible that such strong A-U correlation appears 
only when approaching the phase transition. This can 
be investigated by evaluating the slow A-U correlation 
coefficient as a function of temperature. Also, it is pos- 
sible that such a correlation only appears on even slower 
times-scales (e.g. /is). 



CONCLUSIONS 

This paper reports a study of thermodynamic equilib- 
rium fluctuations of phospholipid membranes. On long 
time scales, we identify strong volume-energy correlations 
of a kind that were previously only observed for simple 
model liquids . These correlations are mainly doc- 
umented for the fluid phase, but we also show that they 
exist in the ordered phase. 

It may seem surprising that a complex system like a 
biomembrane exhibits such strong thermodynamic corre- 
lations. The identification of the origin of the correlations 
as deriving from the van der Waals interactions of the 
hydrophobic part of the membrane, however, points to a 
common origin of strong thermodynamic correlations in 
simple van der Waals liquids and biomembranes. This is 
consistent with the finding that there are strong energy- 
volume correlations in both fluid and ordered phases: 
The correlations do not depend on the degree of chain 
order, just as for simple liquids where the strong correla- 
tions survive crystallization [11] |T2] . We find weaker cor- 
relation between energy/volume and "geometrical" order 
parameters such as membrane area, thickness, and Scd 
(characterizing ordering of acyl-chains). Thus, one pa- 
rameter is not sufficient to describe thermodynamic fluc- 
tuations. 

Regarding the Heimburg- Jackson nerve signal theory, 
our findings largely confirm one crucial assumption of this 
theory, namely that volume and energy (enthalpy) corre- 
late for microstates. We find strong correlations only on 
the nanosecond and longer time scales, which are, how- 
ever, precisely the relevant times for nerve signals. 
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